/* 
 * File:   integracion.h
 * Author: jonk
 *
 * Created on 15 de enero de 2012, 23:43
 */

#ifndef INTEGRACION_H
#define INTEGRACION_H

double integracion(particula *part, parametros param)
{
    double epoten;
    
    epoten= interaccion(part,  param);
    
    for(int i=0;i<param.npart;i++)
    {
        for(int j=1;j<3;j++)
        {
            part[i].pos[j]+=part[i].force[j]*param.dt;
        }
#ifdef COUETTE
        part[i].pos[0]+=(part[i].force[0]+part[i].pos[2]*param.gamma0*param.omega*cos(param.omega*param.t))*param.dt;
#else
        part[i].pos[0]+=part[i].force[0]*param.dt;
#endif
    }
    
    return epoten;
}

#endif  /* INTEGRACION_H */
